Continuum effects and three-nucleon forces in neutron-rich oxygen isotopes 
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We employ interactions from chiral effective field theory and compute binding energies, excited 
states, and radii for isotopes of oxygen with the coupled-cluster method. Our calculation includes 
the effects of three-nucleon forces and of the particle continuum, both of which are important for 
the description of neutron-rich isotopes in the vicinity of the nucleus ^''O. Our main results are the 
placement of the neutron drip-line at ^*0, the assignment of spins, parities and resonance widths 
for several low-lying states of the drip-line nucleus, and an efficient approximation that incorporates 
the effects of three-body interactions. 
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Introduction. - Neutron-rich oxygen isotopes are par- 
ticularly interesting nuclei. First, the nuclei ^^O and ^^O 
exhibit double magicity at the neutron numbers TV = 14 
and N — 16, respectively, see for example Refs. [iHll. 
Second, oxygen is the heaviest element for which the neu- 
tron drip line is established experimentally. The recent 
experiments 0, Q show clearly that the nuclei 25,26q ^^.^ 
unbound, thus making ^'^O the most neutron-rich bound 
isotope of oxygen. The spectroscopy of the drip line nu- 
cleus ^^O was studied in a recent experiment @. One 
of the exciting results of this study is a state with an 
unknown spin and parity at about 7.5 MeV of excita- 
tion energy. Theoretical studies in this region of the nu- 
clear chart are challenging 0-0] ■ Volya and Zelevinsky f?\ 
employed an empirical two-body shell-model interaction 
(above the core of ^^O) and included the particle contin- 
uum in their calculation of neutron-rich oxygen isotopes. 
Otsuka et al. [§| included three-nucleon forces (3NFs) 
within the sd-shell model (keeping ^^O as a core with 
empirical single-particle energies) and found that three- 
body forces yield ^^O at the neutron drip line. The ab 
initio computations of neutron-rich oxygen isotopes by 
Hagen et al. Q employed microscopic interactions from 
chiral effective field theory [l^l, had no core, but were 
limited to nucleon-nucleon (TV TV) interactions. Thus, we 
are still lacking a complete computation of neutron rich 
oxygens that properly accounts for (i) the effects of three- 
nucleon forces, (ii) the presence of open decay channels 
and particle continuum, and (iii) many-nucleon correla- 
tions. It is the purpose of this letter to fill this gap. In 
particular we predict the spins and lifetimes of several 
resonances in ^^O at excitation energies around 7 MeV 
and thereby shed light on the recent experiment [6| . 

Effective field theory (EFT) is the framework that al- 
lows for a consistent formulation of low-energy nuclear 
Hamiltonians and currents jiol - fl^ . Within chiral EFT, 
3NFs are important contributions that enter at next-to- 



next-to-leading order in the power counting. In light 
nuclei, our understanding of 3NFs has improved consid- 
erably over the past decade, and nuclear binding ener- 
gies, spectra and decays cannot be understood without 
them jll-[i3. The study of the role of 3NFs in medium- 
mass nuclei and exotic, neutron-rich nuclei is a frontier in 
contemporary nuclear structure theory. To date, the full 
inclusion of 3NFs is limited to p-shell nuclei. For heavier 
nuclei or nuclear matter, several approaches [^, [isl - Eoj l 
employ a normal-ordered approximation [2l| . resulting 
in a medium-dependent two-body potential that includes 
effects of the 3NFs. Furthermore, the employed interac- 
tions are renormalization group transformations (22| of 
interactions from chiral EFT. 

Hamiltonian and model space. - The intrinsic A- 
nucleon Hamiltonian used in this work reads 



Vnn + V^SN. 



off ■ 
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Here, T is the kinetic energy and Tcm is the kinetic en- 
ergy of the center of mass and depends on the mass num- 
ber A. The potential Vnn is the chiral NN interaction 
developed by Entem and Machleidt fioj at next-to-next- 
to-next-to-leading order (N^LO) within chiral EFT. The 
potential VsNefif is the in-medium NN interaction de- 
rived by Holt et al. (20j from the leading order chiral 
3NF by integrating one nucleon over the Fermi sea (i.e. 
up to the Fermi momentum kp) in symmetric nuclear 
matter. The leading chiral 3NF depends on five low- 
energy constants (LECs). The LECs Ci = -0.81 GeV-\ 
C3 = —3.20 GeV~^, and C4 — 5.40 GeV~^ appear also in 
the two-pion exchange part of the chiral NN interaction 
and have the same values as in the N^LO NN potential 
we employ [10]. The remaining LECs of the 3NF are set 
at CD — —0.2 and ce = 0.7 together with A^^ = 0.7 GeV. 
For the oxygen isotopes considered in this work we ap- 
ply the Fermi momentum kp ~ 1.0 fm^^ in our potential 
VsNcff • Consistent with the NN force, the effective cutoff 
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for the 3NF is A = 500 MeV. 

Let us comment on our two-body potential VaNcff that 
contains the effects of 3NFs. The normal-ordered ap- 
proximation of 3NFs [l^ [in still requires one to com- 
pute an enormous number of three-body matrix elements. 
This poses a great challenge for the large model spaces 
we need to consider. The approach of this Letter is thus 
simpler: The summation over the third particle is per- 
formed in momentum space before the transformation to 
the oscillator basis takes place [38|. This procedure avoids 
the costly computation of three-body matrix elements in 
large oscillator spaces, but it introduces an approxima- 
tion by replacing the mean-field of a finite nucleus by 
that of symmetric nuclear matter. To correct for this 
approximation, we adjusted the LEG ce away from the 
optimal value established in light nuclei (23 |. 

The coupled-cluster method is essentially a similarity 
transformation of the Hamiltonian with respect to a ref- 
erence state. This method is accurate and efficient for 



nuclei with closed (sub-) shells |25l - l27| . We compute the 
ground states of 16,22, 24, 28q -^j^yn the singles and dou- 
bles approximation, while three-particle-three-hole (3p- 
3h) excitations are included in the A-CCSD(T) approach 
of Ref. [2^. For excited states in these closed-shell iso- 
topes we employ the equation-of-motion (EOM) coupled- 
cluster method with singles and doubles. The open- 
shell nuclei i547,21,23,25q g^j-g computed within the par- 
ticle attached/removed EOM formalism, and we employ 
the two-particle attached EOM formalism [s^] for the 
nuclei ^^'^^O. For details about our implementation see 
Ref. (2^ . These EOM methods work very well for states 
with dominant Ip-lh, Ip, Ih, and 2p structure, respec- 
tively. We use a Hartree-Fock basis built in 17 major 
oscillator shells and varied the oscillator spacing huj be- 
tween 24 MeV and 32 MeV. Well converged energy min- 
ima are found at fia; « 28 MeV for all oxygen isotopes. 
Open decay channels and the particle continuum near the 
dripline nucleus ^^O are included within the Gamow shell 
model [3ll [321 ■ The single-particle bound and scattering 



states result from diagonalizing a spherical Woods-Saxon 
Hamiltonian in a discrete momentum basis in the com- 
plex plane [32, 33|. In the case of computing resonances 



in O we used 35 mesh points for the partial wave 
on a rotated/translated contour in the complex momen- 
tum plane as described in Ref. [s^] . The excited states 
we compute in 22,24q g^j.^ dominated by Ip-lh excitations 
and continuum mixing from other partial waves is small. 
In computing radii we discretized the real momentum 
axis with 40 points for the neutron and proton partial 
waves closest to the threshold. This guarentees the cor- 
rect exponential decay of matter and charge densities at 
large distances. 

Results. - Figure [T] shows the ground-state energies 
of the computed oxygen isotopes (red squares) compared 
with experimental data (black circles) and results lim- 
ited to chiral NN interactions only (blue diamonds) . For 



the isotopes around ^^O, NN interactions alone already 
describe separation energies rather well, and the inclu- 
sion of effects of 3NFs mainly changes underbinding into 
overbinding. For the more neutron-rich oxygen isotopes, 
the 3NFs significantly change the systematics of the bind- 
ing energies, and energy differences are particularly well 
reproduced. The nuclei 25,26q ^j-g unbound with respect 
to 240 by about 0.4 MeV and about 0.1 MeV, respec- 
tively, in good agreement with experiments [1[3. We 
predict ^^O to be unbound with respect to ^'^O by about 
4 MeV and with a resonant width of about 1 MeV. The 
extremely short life time of ^^O poses a challenge for ex- 
perimental observation. The energy difference between 
light and heavy oxygen isotopes is not correctly repro- 
duced when compared to data. We believe that this 
is due to the fact that our interaction VsNcff is based 
on symmetric nuclear matter. For smaller values oi kp, 
the ground-state energy of the lighter oxygen isotopes is 
increased (and can be brought to good agreement with 
data), while the heavier isotopes are significantly under- 
bound. The value we chose for kp is thus a compromise. 
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FIG. 1: (Color online) Ground-state energy of the oxygen 
isotope '^O as a function of the mass number A. Black cir- 
cles: experimental data; blue diamonds: results from nucleon- 
nucleon interactions; red squares: results including the effects 
of three-nucleon forces. 

Let us comment on our computation of oxygen iso- 
topes with open shells. First, we solve the CCSD equa- 
tions for the Hamiltonian ([Ij of the closed-shell refer- 
ence state, but employ the mass number A ± 1 in the 
kinetic energy of the center-of-mass. In a second step, 
we add (remove) a neutron within the particle attached 
(removed) EOM. This procedure ensures that the final 
result is obtained for the intrinsic (i.e., translationally 
invariant) Hamiltonian of '^^^O. The J'^ = 1/2+ ground 
state energy of ^^O, shown in Fig. [U resulted from us- 
ing particle-removal EOM from ^'^O. For ^^'^^O, we per- 
formed a two-neutron attached EOM computation based 
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on the reference states for ^^'^^O, the latter being com- 
puted with the mass number A = 18, 26 in the kinetic 
energy of the center of mass. This approach is unprob- 
Icmatic for separation energies but it introduces an error 
in the computation of resonance widths. Our computa- 
tion of ^^O within the neutron attached EOM employs 
a Gamow basis. Here the continuum threshold is incor- 
rectly set by the closed-shell reference of ^''O, computed 
with the mass number A = 25 in the kinetic energy of 
the center of mass. Clearly, this introduces a small error 
by shifting the scattering threshold of the order (Tcm) /A, 
and thereby affects the widths of resonance states that 
are very close to the threshold. In Ref. [s^ we showed 
that the coupled-cluster wave function factorizes into an 
intrinsic and a center-of-mass part. The center-of-mass 
wave function is to a very good approximation a Gaus- 
sian with a frequency hui « 14MeV for ^''O. Thus, we 
estimate the error introduced in the scattering threshold 
of to be about 0.4 MeV. 

Figure [2] shows the excitation spectra of neutron-rich 
oxygen isotopes and compares the results limited to chi- 
ral nucleon-nucleon interactions (blue lines) to the results 
obtained with our inclusion of 3NFs (red lines) , and to ex- 
perimental data (black lines), [s^ The particle continua 
above the scattering thresholds are shown as gray bands. 
Calculations limited to NN interactions yield the correct 
level ordering but very compressed spectra when com- 
pared to data, and all the computed excited states are 
well bound with respect to the neutron emission thresh- 
olds. However, the inclusion of 3NFs increases the level 
spacing and significantly improves the agreement with 
experiment. Several of the excited states are resonances 
in the continuum, and the proximity of the continuum is 
particularly relevant for the dripline nucleus ^^O. Here, 
the Gamow basis is essential for a proper description of 
the excited resonant states. In ^^O we find three resonant 
states near the unknown experimental state at about 
7.5 MeV of excitation energy [1]. The excited J"" = 3/2+ 
state in ^'^O is computed as a neutron attached to ^^O, 
while the excited = 5/2+state is computed by neutron 
removal from ^^O. Since we are interested in the excita- 
tion energy relative to the ground state we compute the 
J'^ = 1/2+ ground state either by adding or removing a 
particle, consistent with the particular excited state. For 
the lighter isotopes i^-i^'i'^O, our inclusion of 3NFs yields 
only smaller changes to the spectra when compared with 
results from NN interactions only. 

For the closed-shell nucleus ^"'O we also computed res- 
onance widths of excited states that are dominated by 
Ip-lh excitations from the Si/2, ^5/2 hole to the ^3/2 par- 
ticle orbitals. Table U shows that the first J'^ = 2+ 
and 1+ excited states agree well with experimental data, 
both for the excitation energy and the resonance widths. 
The rather small widths and quasi-bound nature of these 
states can be attributed to the large angular momentum 
barrier of the ^3/2 orbital, together with neutron pairing 
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FIG. 2: (Color online) Excitation spectra of oxygen isotopes 
computed from chiral nucleon-nucleon interactions, with in- 
clusion of the effects of three-nucleon forces, and compared to 
data. 
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TABLE I: Excited states in ^"O computed within EOM- 
CCSD compared to experimental data from Ref. 01 . Energies 
and widths are in MeV. 



effects. Above these 2+ and 1+ resonances we find sev- 
eral states with with spin and parity J'^ = 1+ to 4+ and 
excitation energies ranging from 6.8 MeV to 8.9 MeV. 
The small ratio /E^+ w 1.31 and the relatively high 
energy lend theoretical support to the doubly magic 
nature of ^^O The low experimental resolution of the 
resonance at 7.5 MeV let Hoffman et al. [ij to speculate 
that this resonance could be a superposition of narrow 
resonances with spins and parity J'^ = 1+ to 4+. Our 
calculation clearly supports this suggestion, except for 
the I2 state which we find at 8.9 MeV of excitation en- 
ergy and with a width of 0.57 MeV. 

Figure [3] shows the computed point matter and point 
charge radii for the neutron-rich isotopes 2i-24q -^^i]^ 
comparison to the experimental data (Ref. [s^ for -^^O 
and Ref. IS?] for ^^-^^O). This computation employs the 
intrinsic density with respect to the center of mass. Our 
radii agree fairly well with experiment, and for 22,23q 
they are within the estimated error bars. The com- 
puted charge radii clearly exhibit an odd-even stagger- 
ing consistent with the shell closures at neutron num- 
bers N = 14, 16. For ^^O, our computed point matter 
and charge radii are 2.27 fm and 2.29 fm, respectively. 



4 



Compared to experiment our radii are about 0.25 fm too 
small, this is due to the overbinding and too large neu- 
tron and proton separation energies we compute for ^^O. 
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FIG. 3: (Color online)Point matter and point charge radii of 
neutron rich oxygen isotopes computed from chiral nucleon- 
nucleon interactions, with inclusion of the effects of three- 
nucleon forces, and compared to data. 

Summary. - We employed interactions from chiral ef- 
fective field theory, performed coupled-cluster computa- 
tions of oxygen isotopes, and included effects of the parti- 
cle continuum and of three-nucleon forces. Three-nucleon 
forces were approximated as density-dependent nucleon- 
nucleon forces. This approach is computationally feasi- 
ble and in keeping with the spirit of effective field theory. 
Compared to computations based on nucleon-nucleon in- 
teractions alone, the included 3NFs yield a significant 
improvement in binding energies and spectra. Our re- 
sults confirm that chiral interactions yield the neutron 
drip line at ^^O, and we are able to compute spin, pari- 
ties and resonance widths for several excited states close 
to the dripline. In particular, we compute several long- 
lived resonances at about 7 MeV of excitation energy in 

240. 
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